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ABSTRACT 

We present the analysis of an IRS 5-38 /im spectrum and MIPS photometric measurements of an 
infrared echo near the Cassiopeia A supernova remnant observed with the Spitzer Space Telescope. 
We have modeled the recorded echo accounting for PAHs, quantum-heated carbon and silicate grains, 
as well as thermal carbon and silicate particles. Using the fact that optical light echo spectroscopy has 
established that Cas A originated from a type lib supernova explosion showing an optical spectrum 
remarkably similar to the prototypical type lib SN 1993J, we use the latter to construct template 
data input for our simulations. We are then able to reproduce the recorded infrared echo spectrum 
by combining the emission of dust heated by the UV burst produced at the shock breakout after the 
core-collapse and dust heated by optical light emitted near the visual maximum of the supernova 
light curve, where the UV burst and optical light curve characteristics are based on SN 1993J. We 
find a mean density of ~680 H cm~'^ for the echo region, with a size of a few light years across. We 
also find evidence of dust processing in the form of a lack of small PAHs with less than ^300 carbon 
atoms, consistent with a scenario of PAHs destruction by the UV burst via photodissociation at the 
estimated distance of the echo region from Cas A. Furthermore, our simulations suggest that the weak 
11 /im features of our recorded infrared echo spectrum are consistent with a strong dehydrogenated 
state of the PAHs. This exploratory study highlights the potential of investigating dust processing in 
the interstellar medium through infrared echoes. 

Subject headings: ISM: clouds - ISM: lines and bands - ISM: evolution - supernovae: individual (Cas 
A) - Infrared: ISM 



1. INTRODUCTION 

Light echoes refer to light emitted by a bright source, 
such as a supernova (SN), which is scattered by dust 
located off the direct line-of-sight. The additional path 
length from the SN to the dust cloud and then to the 
Earth induces the time delay of the echo. Depending on 
the system considered and the brightness of the source, 
light echoes can be seen several hundreds of light years 
away from the source, and subsequently several hundred 
years after the initial SN outburst observation. 

Infrared (IR) echoes differ from optical light echoes in 
that the SN light is not scattered, but rather absorbed 
by dust, which will the n re-radiate it at IR wavelengths. 
iBode and EvansI (J1980I ) considered this effect as the most 
plausible explanation for the IR excess see n in SN 1979c , 
a point which was later on confirmed bv iDwed (|1983f ). 
IR echoes spanning over an area of several square-degrees 
on-sky have been discovered around the youngest known 
gal actic core-c oUapsc SN remnant Cassiopeia A (Cas 
A) (jKrause et al. 2005 ). Dwek and Arendt (2008) used 
their spectra to deduce that these must have b een created 
by a UV flash. The same echoes also enabled iKim et all 
|2008) to reconstruct a 3D map of the ISM surrounding 
Cas A, a concept that was previously used with SN 1987A 
by Xu and Crotts (1999). 

The Cas A supernova remnant is believed to be the 
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remains of a stripped-envelope type lib core-collapse 
from a red supergiant progenitor with an initial mass 
in the range of 15-25 Mq which might have lost much 
of its hydrogen envelo pe due to a bi nary interaction 
(|Young et al. 2006 : Krauseet al.ll2008D . and is located 
at R.A.: 23''23'°27.77" ; Dec: -h58°48'49.4" (J2000) 
(jThorstensen et all 1200 1[ ). Studies of the proper motion 
of ejecta with the Hubble Space Tele scoye suggest tha t 
the explosion occurred in 1681±19 (Fes en et al.l 120061. 
Sir John Flamsteed (1646-1719) might have visually ob- 
served the supernova on 1680 August 16, but it is not 
clear whether his myster ious 3 Cassiopeia observation 
was actually Cas A or not (jAshwor th 1980; Hughes 1980). 

Along with the infrared echoes, several optica l light 
echoe s have b een discove r ed aro un d Cas A (Krau s e et al] 
l2005l I200I IRest et all [200l . IKrause et al] (|2008l ) 
showed that the spectrum of one of them is strikingly 
similar to the spectrum of SN 1993J, which enabled the 
precise spectroscopic classification that Cas A was a type 
lib SN. The knowledge of the exact type of the Cas A 
supernova will play a central role in our analysis. 

SN 1993J exploded in M81 in late March 1993. 
It has been extremely well monitored across 
the entire spectrum from its very first moments 
(iFilippenko and MathesonI 120031 ). including in the UV 
(e.g. iJefferv et alj|1994[ ). the wavelength range at which 
the irradiation spectrum has critical consequences on 
interstellar dust. Such an amount of data makes it a 
target of choice for supernova simulations, that are able 
to simu late the observed spe ctrum already from its early 
phases (jWooslev et al.lll994l ). SN 1993 J experienced an 
initial, short-lived (few hours) intense UV burst, later 
on followed by a longer (several days) less intense optical 
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component (|Blinnikov et al.lll998D . Hence, it is a very 
fortunate situation for our study that Cas A was of the 
same type and similar to SN 1993J - one of the brightest 
SNe in the northern sky for the last century, and for 
which a nearly unparalleled amount of data (including 
UV) exists. 

Infrared echoes can be used to study the structure and 
compositi on of the the rmally radiating dust, as was sug- 
gested by ID wed (|1985D for circumstellar dust shells. In- 
frared echoes especially allow to study several dust pro- 
cessing mechanisms, such as photodissociation, dehydro- 
genation and ionization, in-situ. Due to the fact that 
more than 320 years have passed since Cas A exploded, 
the dust giving rise to the infrared echoes observed today 
is at least 160 light years away from Cas A. These echoes 
therefore most likely originate from interstellar rather 
than circumstellar dust. Such an analysis of infrared 
echoes is similar to what ISugermanI (20031 suggested for 
optical light echoes, but using their infrared counterpart. 
As the SN light is absorbed and re-radiated at longer 
wavelengths, the resulting IR spectrum is a complex sig- 
nature of both the SN light and the dust composition and 
characteristics. Using optical light echoes to identify the 
SN type, one can create a template light curve. Feeding 
this template in a suitable dust modelling program, one 
can simulate the IR echo spectrum - and use the result- 
ing simulations to understand how ISM dust is behaving 
when illuminated by the SN light. 

One of the main advantages of this method is that IR 
light echoes potentially give us access to the consequences 
of pristine interstellar material subject to a strong and 
short radiation burst. As such, assuming an undisturbed 
initial ISM, the dust signature detected in the IR light 
echo will be the consequence of the SN burst only, as no 
other mechanism, such as dust formation for example, 
can be acting on the very short timescale of a SN burst. 

In this article, we explore the potential of this ap- 
proach by analyzing an IR echo spectrum recorded with 
the Spitzer Space Telescope around Cas A. We describe 
our observations and the resulting recorded IR echo spec- 
trum in Sec. [2J Wc discuss how we simulate a light echo 
in Sec. [31 introduce our dust model in Sect. 21 and present 
our Cas A SN template burst spectra in Sec. O Our re- 
sults are detailed in Sec. [HI In Sec. [71 we discuss how 
PAHs destruction (Sec. 17. 2|) . ionization (Sec. 17. Sp and 
dehydrogenation (Sec. 17.41) can improve our fit to the 
recorded IR echo spectrum. We comp are our results with 
the work of IDwek and Arendtl (|2008D in Sec. M Finally, 
we explore the validity of our dust destruction hypothesis 
in Sec. [HI and present our conclusions in Sec. [TUl 

2. OBSERVATIONS AND DATA REDUCTION 

The observations used in this study were obtained 
with instruments ab oard the Spitzer Space Telescope 
(jWerner et al.l I2004D . After the first discover y of in - 
frared echoes around Cas A by iHines et"al] (|2004D . 
regular monitoring observations were carried out with 
the Multiba nd Imaging Photometer for Spitzer (MIPS; 
IRieke et al.l (|2004)) (PIDs: 233, 20381, 30571; Pis: 
G.Rieke, O.Krause). A roughly 3x3 square degrees 
field was observed every year and regions with especially 
interesting activity were observed at an interval of about 
half a year. MIPS observations were re duced using t he 
MIPS Data Analysis Tool (DAT; .Gordon et al.l ([2005I) ). 



starting at the RAW data product level. 

Using those monitoring observations, a particularly 
bright source was detected at the location R.A.: 
23'^21™40"; Dec: -h59° 34' 25". Fig. [T] shows its evo- 
lution at 24//m over a time period of 3 years. It lies at 
an angular distance of 2857 arcsec away from the optical 
expansion center of Cas A. Although the source is already 
visible in the first epoch of observations, it changes its 
morphology, brightness and position drastically within 
just a few months time. In early 2006, it appears that 
most emission comes from three point like sources, sug- 
gesting that the real echoing structure is still unresolved 
with Spitzer. The emission diminishes very quickly be- 
tween February and October 2006 indicating that the SN 
light beam swept past a probably very well defined end 
of the echo region. 

On 2006 January 29, a low- resolution spectrum of this 
bright echoing sou rce was obta. i ned w ith the Infrared 
Spectrograph (IRS: iHouck et all ()2004[) ) as part of pro- 
gram (PID) 20381 (P.I.: O. Krause). A 24 fim aerial 
image of the echoing source is shown in Fig. [2l where 
we have overlaid the IRS red peak-up insert (rectangle) 
as well as the spectroscopic slits position; 1 x Long-Low 
(LL) and 3x Short-Low (SL). The red peak- up function- 
ality of IRS was used to achieve the best possible center- 
ing of the slits on this fast moving echo-feature. In ad- 
dition to space-based infrared observations, we also ob- 
tained ground based near-infrared (A=2.2/im) Ks-band 
images for this echo region using Omega 2000 on the 
Calar Alto (Spain) 3.5m telescope in December 2005, 
which is shown in the bottom of Fig. [5] for comparison. 
These images show the similar overall structure of the re- 
gion, but reveal very small scale filamentary substructure 
of the echo region. A full analysis of these observations is 
to be presented in Besel et al. (in preparation) . It should 
be noted here that the Ks-band as well as the J-band and 
H-band images (which have similar morphology, and are 
not shown here) show the scattered light echo, not the 
IR radiating dust echo. 

The S18.7.0 version of the Basic Calibrated Data 
(BCD) products from the Spitzer Science Center (SSC) 
pipeline were fed into SMART (jHigdonet al.l l2004t 
iLebouteiller et alll20100 for final reduction. Depending 
on the used module, total exposure times ranged between 
120s and 300s. One spectra was extracted for each of the 
two brightest point-like sources. They were averaged sep- 
arately and combined later. A background was estimated 
based on the non-echo-contaminated parts of either the 
off-source order or the individual slits (contamination is 
easily evident on the MIPS 24 /zm image for one end of 
the LL slit). As the SL and LL modules cover slightly 
different special regions of the echoing source the result- 
ing spectra were scaled to match each other's continuum 
levels. Absolute calibration was verified by comparing 
the integrated IRS flux and the measured 24 /im pho- 
tometry. 

Fig. [3] shows the final combined extracted spectrum of 
this IR echo. Also plotted are MIPS 24 fim and 70 iim 
photometric measurements within a 22.2 arcsec aperture 
(white squares). Note that their error bars are of the or- 
der of the symbol size. Black rombuses at 5.7 /xm, 7.9 /Ltm 
and 24 /j,m correspond to the Spitzer Synthetic Photo- 
metric points of the echo spectrum, reconstructed using 
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Fig. 1. — Time evolution of the echo region analyzed in this study. Observation date and Spitzer Program Identification (PID) are 
imprinted in each panel. All observations are MIPS 24^m images, except for an IRS red-peakup image (18.5-26 ^m) on 2006 January 29. 
Each panel is displayed in logarithmic scale, with North up and East towards the left with a field of view of roughly 51" X 69". 
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Fig. 2. — Top: 24 fim aerial image of the echoing source with IRS 
red-peakup insert (dashed rectangle) and footprints of the Long- 
Low (LL, in green in the online version) and Short-Low (SL, in 
blue in the online version) spectroscopic slits positions. Bottom: 
Calar Alto Omega 2000 Ks-Band image of the same area. Visible 
are tiny structures within the infrared emitting area. 

the Recipe 10 of the Spitzer Data Analysis Cookbook 0. 
As the IRS slits do not cover the entire echo region, the 
spectrum and synthetic photometric points were scaled 
to match the 24 /im photometry of the entire region. 

3. SIMULATING THE INFRARED ECHO 

The spectrum and photometric points presented in 
Fig. [3] are the result of the interaction of the Cas A SN 
burst and a region of the ISM located away from the 
SN center. In that respect, the data potentially contain 
information on both the SN explosion and the echoing 
dust. A simple visual inspection of the spectrum already 
reveals a very weak 11 /im complex of PAH features as 
compared to the strength of the 7 /im PAH features, a 
point that will be discussed further in Sec. 17.41 

Here, we adopt the following strategy to study and ex- 
tract the underlaying information on the echoing dust 
contained in the measured echo: we implement a dust 

* see |http://ssc. spitzer. caltech.edu/dataanalysistools/cookbook/ 1 
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Fig. 3. — Final extracted IRS spectrum of the echo region (in red 
in the online version) and associated Spitzer synthetic photometric 
measurements (black rombuses). The white squares correspond to 
MIPS 24 /xm and 70 /^m photometric measurements of the same 
echo region. The spectrum and synthetic photometric points have 
been scaled to match the 24 ^m MIPS photometry of the entire 
region. 

modelling program (see Sect. |4]) to simulate the IR re- 
sponse of ISM dust when heated by radiation with a cer- 
tain spectral energy distribution (SED). To obtain the 
dust thermal emission response, its distribution along the 
line of sight has to be convolved with the corresponding 
irradiation spectrum from the supernova. We assume 
that the effective radiation field can be decomposed by 
(I) the UV dominated SED during the first hours after 
the shock outbreak and (II) the optically dominated SED 
at the epoch near maximum light (see Sect.[5|). We adopt 
those two template SEDs to account for the convolution 
of the input radiation spectrum with the echo region on 
length scales of at least a few light days. 

In Fig. m we present a geometrical sketch of the echo 
principle. The position of the echo region is located, at 
any given ti me, on a ellips oid with Cas A and the Earth 
at the foci ()Couderdll939() . The distance r between the 
echo region and the center of Cas A is : 



d{d + ct) (l - cos(a)) + {ctf 12 



d{\ 



cos 



(a)) ^' ct 



(1) 



where t is the time delay between the light of Cas A's 
explosion reaching Earth and the observation of the light 
echo, d the distance to Cas A and a the on-sky angular 
separation between Cas A and t he echo region. U sing 
a distance to Cas A of 3.4 kpc (iReed et alJ Il995f) and 
an explosion date of AD 1681 (jFesen et alJ 120061 ) . the 
absolute distance of the echo region to the supernova 
remnant is r = 199 light years. 

A schematic of the echo region itself is shown in Fig. [5j 
Because the UV burst and optical phase do not occur 
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Fig. 4. — Light-Echo geometry. The observer on Earth and the 
Supernova are located in the foci of an eUipsoid. All points on 
this ellipsoid describe locations with equal delay time t which cor- 
responds to the time difference between the original SN outburst 
observation and the current observing epoch. Material situated on 
any point of the ellipsoid can simultaneously be seen as an echoing 
source. 

simultaneously, they will heat a given area of the echo 
region sequentially. In Fig. [5l area (A) corresponds to 
the echo region not seen excited at the moment of the 
observation, while areas (B) and (C) denote respectively 
the schematic location of the UV burst and optical light 
curve excited slabs, moving away from Cas A. 
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Fig. 5. — Schematic of the echo region. The echo emission is 
coming from two defined regions (B) and (C) heated by the optical 
and UV burst respectively. The drawing depicts what we see at 
a given time from the Earth, which is different from the state 
of the echo region at a given time. With a thickness of a few 
light days, the excited slabs are much smaller than the actual echo 
region, which is not completely excited at one given moment in 
time (region (A)). 

It is important to note here that if we see those ex- 
cited slabs simultaneously, they are not excited simulta- 
neously, because of the finite light speed. The UV-burst 
excited slab (a.k.a. UV-heated slab) is located farther 
away from Cas A and the Earth, and must therefore 
be excited earlier than the optically-excited slab (a.k.a 
optically-heated slab) to be seen at the same time on 
Earth. In situ, the two SEDs will be seen by a dust 



grain with the same time interval as when they occurred 
after the shock breakout. As seen from the Earth, the 
excited slabs are encompassed in between ellipsoids of 
different radii rg, ri, r2 and r^ defined by Eq. [1] with io, 
ii , ^2 and ^3 the delay times between the observation and 
the shock breakout, the end of the UV burst, the start 
of the optical phase and the end of the optical phase, 
respectively. 

Let us define isky,|| and isky.i. the on-sky dimensions 
of the echo, in the direction of and perpendicular to Cas 
A, respectively. Because the UV-heated and optically- 
heated slabs thickness will only be of a few light days, 
the excited volume of the echo region visible at one given 
moment in time is much smaller than the total volume 
of the echo region. For simplification, because we do not 
resolve individual ISM structures and because the echo 
region is small compared to the ellipsoid and its local cur- 
vature depicted in Fig. 21 let us assume the excited slabs 
to be parallelepipeds, of which the sections are depicted 
by the light-grey (B) and dark-grey (C) regions in Fig. [5] 
(in orange and blue in the online version). Their height 
(perpendicular to the drawing plane) is Lsky._L- Their 
respective thickness Ljjy and Lopt along the line-of-sight 
are : 



where, from Fig. 2] and the law of cosines 



-^uv = (a;o -2^1) (2) 

Lopt = {x2 - X3) (3) 

he law of cosines : 

Xi^ct. + d-n Vi = 1,2,3,4 (4) 

For the echo considered here and presented in Sec. [U 
xo = 3439 pc. We can now define an upper limit to the 
visible volume of the echo region, T4cho,max, as seen from 
the Earth at any given time : 

K;cho,max — -^sky,|| X (^UV + -^opt) X -^sky,± (5) 

where isky.|| x Lopt and -Z^sj^y || x Lyv are the respective 
areas of the regions (B) and (C) in Fig. \5\ The shape 
of a fictional echo region (dotted surface) depicted in 
Fig. [5] illustrate why 14;cho,max is an upper limit on the 
visible volume of the echo region, as the dust does not, 
in principle, have to fill completely the areas (B) and 
(C). We delay the numerical computation of Vocho,max 
until Sec. 16.11 after introducing our adopted template 
UV burst and optical light curve in Sec. [SJ 

4. SIMULATING ISM DUST 

Simulations provide an efficient tool to study the the- 
oretical behaviour of interstellar dust in an attempt 
to unveil its composition and detailed characteristics. 
Various models^ such as the silicate core-carbonaceous 



("Des ert et al.lll990l : IJones et a l. 1990; Li and Grccnberg 
1997p, the composit e (IMathis' and Whifl:cn.l989.:.Mathi£ 
19961 iZubko et al] I2004D or the silicate-graphite- 
Polycyclic Aromatic Hydrocarbon s (PAHs) model 
Draine and Lcc 1984; Sicbcnmorgcn and Kriigel 119921 : 



Li and Drainell2001l : iDraine and Lill2007[ ) have been pro- 
posed to reproduce the behaviour of ISM dust. To study 
the recorded IR echo spectrum shown in Fig. [31 we have 
implemented our own mod elisation progr a m, st rongly 
based on the dust model of iLi and Draind ()200lD . Our 
dust mix contains PAHs, quantum-heated carbon and 
silicate grains, as well as thermal graphite and silicate 
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grains. Similarly to lLi and Draind (|2001[) . the PAHs and 
quantum-heated carbon grains are mixed and form what 
is re ferred to as the quaii tum- heated carbonaceous dust. 
The lLi and Draing ()2001t ) model is a well established and 
observationally tested model for the average dust com- 
position of the general Milky Way interstellar medium. 
Since the echo region considered in this study is located 
in the general ISM at a distance of ~199 light years from 
the supernova, using this model appears very appropri- 
ate. 

We simulate particles from ^3.5 A up to ^6000 A, 
and have a transition from quantum-heated to thermal 
grains a t ~250 A. The dust mix is fol lowing the for- 
mula of iWeingartner and Draind (I2001D m ultiplied by 
0.92, as suggested bv lDraine and Lil (|2007[ ). from whom 
we have implemented the updated Drude profiles for the 
PAHs opacitie s . Ou r simulation yet differs from the 
iLi and Draind (I2001D model in a few aspects^ We use 
alSieber 



enthalpies from lSiebenmorgen and K riigol (1992Dj which 
are based on a formula bv IChase et al.i (|1985i ). Our sil- 
icate opacities are calculated via the Mie t heory using 
scattering and absorption cross section from llvezic et al.l 
([1997) . Finally, we also use enthalpies of graphite for the 
silicate grains for simplification. 

Quantum-heated grains are subject to stochastic heat- 
ing by the incoming photons. They do not behave as 
constant-temperature grains, but experience strong tem- 
perature variations, and a careful approach is required 
to compute their emission. More precisely, their emis- 
sivity is obtained by calculating the temperature prob- 
ability distribution for a certain grain size, and we re - 
fer the reader to t h e work of [Praine and Lil (J2001I ): 
ILi and Drai'^ (|2001D : IDraine and Lil (J2007D fo7 detailed 
explanations regarding the computation of the emissiv- 
ity of small grains. Large thermal grains, on the other 
hand, will experience much smaller temperature excur- 
sions, and can be approximated as being in thermal equi- 
librium. Under this assumption, their emissivity is ob- 
tained by comparing the incoming radiation to the re- 
radiated energy dLi and Drai nc 2001) : 



abs(a,A)FA(A)dA= / CMa,X)47:Bx{T)dX 



(6) 

with -Fa (A) the incoming radiation flux at wavelength A, 
T the equilibrium mean temperature of the grain, Bx (T) 
the Planck function and Cabs (a, A) the absorption cross- 
section of a grain of radius a at the wavelength A. 

In Fig. ini we present our simulation (full line, in green 
in the online version) of the ISM dust IR emissivity 
when heat ed by the stan dard Interstellar radiation field 
(ISRF) of iMathis et al.l ^Mj (see Fig. El- The con- 
tribution of each individual dust components is plotted 
in dashed (quantum-heated grains) and dotted (thermal 
grains) lines; (a) - quantum-heated carbonaceous dust, 
(b) - quantum- heated silicate dust, (c) - thermal graphite 
grains, (d) - thermal silicates grains. The simulation code 
returns the emissivity of the dust mix per H nucleon, j^, 
in units of erg s~^ sr~^ H~^ Hz~^, and in Fig.[6l we plot 
the spectrum in the form of fj^, in units of erg s~^ sr~^ 
H-i. 

We validate our code output against the results of 
the JLJ and Drainc (200 1.) code, and the comparison is 




A [/U,rn] 

Fig. 6. — IR emissivity of interste llar dust per H nucleon heated 
by the ISRF of lMathis et al.l l|1983l ). assuming a silicate-graphite- 
PAHs composition. Dashed and dotted curves show the individ- 
ual contribution of the quantum- heated carbonaceous dust (a), 
quantum-heated silicate dust (b), thermal graphite (c) and thermal 
silicate grains (d). 

made in Fig. [71 Our HGL simulation (full line, in green 
in the online v ersion ) reproduces closely the one from 
ILi and Draind (|2001| ) (dotted line), with the exception 
of the PAH features, t hat have been update d to repro- 
duce the simula tions of [D rainc and Lij ([2007D . Note that 
the spectrum of lLi and~Draine (200 li) has been scaled by 
0.92 to allow for comparison. 




Fig. 7. — IR emissivity of interstellar dust per H nucleon for 
HGL dust (full line, in green in the online version) and for ISM 
dust illuminated by 10* times the ISRF (dashed line, in blue in 
the online version). The latte r has been scal e d dow n by 10^. The 
corresponding simulations of ILi and Draind l|20011 'l are shown in 
dotted lines. 

The strategy to simulate the IR echo spectrum pre- 
sented in Sec.[5|is to use our cod e with a temp l ate in put 
spectrum, replacing the ISRF of lMathis et al.l (|1983| ) by 
a theoretical SN burst spectrum. In order to validate 
that our code is stable with a more intense heating radi- 
ation field, we simulate ISM dust heated by a radiation 
field equal to 10^ times the ISRF, and compare o ur re - 
sult with the similar simulation of ILi and Draind (|2001D 
(see Fig. 13 in their article). The comparison is shown 
in Fig. |7| (dashed line, in blue in the onli ne version). It 
is clea r that we reproduce very closely the lLi and Draing 
(|2001f ) simulations for stronger radiation fields, with ex- 
ception to, once again, the intrinsic differences due to the 
update of the Drude profiles for the PAHs, which match 
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the simulations of iDraine and Lil (|2007[ ). We conclude 
that our simulation code is stable to intensity variations 
of the incoming radiation field. 

5. MODELLING THE CAS A SN OUTBURST 

Simulating the spectrum of an infrared echo requires 
the knowledge of the exact radiation SED heating the 
dust. Generally, this is not known with precision, as the 
SN light swept past Earth many years before a light echo 
observation. Fortunately, as we mentioned in Sect. [H 
SN 1993J has been found to be a very good template 
for Gas A, and has been the subject of many observa- 
tions and simulations over a large range of wavelengths. 
Hence, following the characteristics of SN 1993 J, we make 
the following simplification mentioned previously: we 
consider the Gas A SN explosion to be made of two dis- 
tinct phases - a UV burst and an optical light phase. 
Each phase consists of a single SED, which we take to be 
representative of the actual spectrum variations over the 
phase duration. 

To create th e template UV burst s pectrum, we use the 
simulations of lBlinnikov et aLl (|1998f ) . In this model, the 
UV burst is close from its maximum luminosity during a 
few hours. The total UV burst lasts a few days, with its 
luminosity reaching a minimum after about ^^10 days. 
Energetic UV photons have a strong effect on dust and 
especially on PAHs, as we will discuss later on. It is 
therefore critical to have a representative SED that con- 
tains those energetic UV photons. The peak of the UV 
burst is however fairly brief, and one has to pay attention 
not to over-estimate the amount of very energetic pho- 
tons. With these considerations in mind, we adopt our 
template Gas A UV-burst to be a black-body spectrum, 
with a temperature of Tuv = 1-5 x 10^ K. In an attempt 
to provide a slightly more realistic light curve, we subse- 
quently adapt the s hape of the black - body curve to that 
of the '.278 curve of iBhnnikov et al.l (|1998[ ) in the 0.001 
- 1 /xm range. We assume a burst duration of 2.5 days, 
and a total luminosity of Luv-b = 6.4 x 10''^ erg s~^. 
For comparison, the luminosity of the black body curve 
prior to the shape modification is Luv=5.9xl0'*'* erg s~^. 
By modifying the burst shape rather than increasing the 
burst temperature, we ensure that we do not overesti- 
mate the number of very energetic photons, while still 
taking into account the fact that a significant amount of 
energy is radiated away during a short period of time. To 
account for the smaller, but yet non- negligible, burst lu- 
minosity beyond the 2-days mark, we adopt a total burst 
duration of 2.5 days, shorter than the ~ 10 days mark, 
to avoid an obvious over-estimate. 

We concede that there exist many other ways to obtain 
a Gas A template UV burst from the data of SN 1993J, 
but any solution, including our adopted template, would 
involve quite a number of assumptions. We emphasize 
here that we do not claim to have implemented neither 
the best, nor the most accurate reproduction of Gas A's 
UV-burst. But rather made an good educated guess. In 
Section [HI we will co mpare our UV b urst template with 
previous estimates bv lDwek and Arendt (2008 ) and show 
that both UV burst esti mates are consist e nt. 

We use the work of [Richmond et al.l (|1994[ ) to cre- 
ate the template optical light curve. They find the 
optical phase to last some 20-30 days, with a lumi- 
nosity depending on the color excess E(B-V) on the 



line of sight. They provide two estimates, a low one 
with E(B-V) = 0.08, and a moderate one, with E(B - 
V)=0.32. Based on the work of iGlocchiatti et al.l ()1995D . 
that finds E(B-V)=0.25 ± 0.05, and using the discussion 
on the absorption on the line-of-sight towards SN 1993J 
of Matheson et al. (2000) , that adopt in their work a 
value of E(B-V) = 0.19, we assume our optical light 
curve to have the characteristics of the [Richmond et al.l 
()1994[) estimates (see Table 15 in their paper) adapted 
for a color excess E(B-V)=0.20. Our optical light curve 
is a black-body spectrum, with a temperature Topt — 
9.4 X 10^ K (between the ^ 7 . 9 x 1(3 ^ and - 11.0 x 10^ 
estimates of [Richmond et all ([1994[ )). a luminosity of 
L^pt ^ 3.8 X lO'*^ erg s'^ ( between the ^ 1 7 x 10 '^^ and 
~ 7.1 X 10"*^ estimates of [Richmond et al.[ ([1994[)) and 
a dur ation of 25 days fsee Fig. 10 in [Richmond et al.l 
(|1994[ )). We stress again that we consider this optical 
light curve as a good educated guess rather than the 
perfect estimate. 

The resulting UV burst and optical light curve fiuxes 
I'Fi, shining on the dust in units of erg s^^ cm~^ are 
shown in Fig. (S) The final UV burst spectrum used 
in our simulations, of whic h the shape has been up- 
dated based on the work of [Blinnikov et al.[ (|1998() , is 
labelled 'UV-B' (dashed line), while the black body UV 
burst with T=1.5xl0^ K is labelled 'UV (dot-dashed 
line) and is plotted for comparison. Using the IDlJj ism- 
tau.prqj routine, we have taken into account some fore- 
ground absorption for the UV burst assuming an atomic 
H-column density Nh of 10^* cm~^, an H II / H I ra- 
tio of 0.26, an He I - and He II - column densities of 
0.1 Nh, and used the ph o toioni zation cross-section of 
[Morrison and McGammonI ([1983[ ). Note that we have 
not taken into account any intrinsic extinction by the 
echo region itself. The various burst characteristics are 
summarized in Table [1] 




Fig. 8. — The UV-B (dashed hne, in blue in the online ver- 
sion) and optical (dotted line, in red in the online version) light 
curve used to simulate the recorded IR echo spectrum shown in 
Fig. \3\ The UV light curve (dot-dashed) is a black body with 
T=1.5xl0^ K and is shown for comparis on, along with the ISRF 
curve (full line) from lMathis et al.l 0M ). 

^ Interactiv e Data Language 

^ see, e.g. http://hea-www.harv ard.edu/PINTofALE/pro/ismtau.prol 
(accessed on 'ZiML December li) 
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6. RESULTS 

6.1. Fitting the echo spectrum 

With the adopted UV burst and optical light curve 
characteristics defined in the previous section, we can 
now get the estimates for the size of the echo region. 
Specifically, with to = 325 year, to ~ ti — 2.5 days, ti — 
i2 = 10 days, t2 — t^ — 25 days, a — 2857 arcsec and 
^sky,|| = -^sky,_L — 1-2 light ycars (22" on-sky aperture, 
see Sec. [U, we get, using Eq. [2]to Eq. \5\ : 

Luv — 1-53 light days 
Lopt = 15.3 light days 
V;cho,max = 6.4 X 10^2 (light ycars)^ (7) 

^5.5 X 10^2 cm^ 

Because we do not resolve the echo region, we make the 
assumption that Fccho the recorded echo flux is equal to 
the sum of the UV-heated slab flux and the optically- 
heated slab flux Fuv and Fopt respectively : 

-Focho = FlJV + Fopt (8) 

Furthermore, Fuv and Fopt relate to the UV-heated and 
optically-heated dust emissivities j,y,uv and jt/,opt in the 
following way : 



Fvv = Ji/,uv X 

-^opt — Ji^,opt X 



Afdust,UV 
?7lH,dust • x'^ 

-Mdust,opt 
"iH.dust • X'^ 



(9) 
(10) 



where Mdust,uv and A/dust, opt are the echoing dust mass, 
in g, heated by the UV burst and the optical light curve 
respectively, mn.dust is the dust mass per H nucleon in 
g H~^, and x is the distance from the Earth to the echo 
region in cm. We ass ume r7iH,dust — 1-89 x 10~^^ g H^^ 
(|Li an d Draind [2001t ). From our assumptions in Sec. |3] 
(see also Fig. [5]), we can express Mdust,uv and Afdust,opt 
as a function of the cloud density Pcioud, and the volume 
of the excited slabs (see Eq. [5] and ^ : 

Mdust,UV = /Ooclio X Lsky,i X isky,|| X {xq ~ Xl) (11) 
A^dust.opt = Pocho X 5sky,± X isky,|| X {x2 - X3) (12) 

Hence, using Eq.[8]to[T2l it is possible to express the echo 
spectrum i<"ocho as a function of the total mass of echoing 
dust Mdust = Af^dust,uv + A^dust,opt in the following way 



-Fccho = [C • >,uv -I- (1 - C) • >,o: 



M, 



dust 



ptj 



"lH,dust • X^ 



(13) 



where Mdust is in g and C = 



A/d 



A/d 



is the ratio of the 



mass of the UV-slab dust to the total echoing dust mass, 
which, using Eq. [TT] and [T^l can be written as : 



c 



Xo -Xl 



0.09 



(14) 



{xo - Xl) + {X2 - X3) 

Eq. [m and [H reflect the fact that the UV-heated 
and optically-heated dust emissivities are mixed with a 
9%/91% ratio (similar to the UV vs optical phase dura- 
tion ratio) . For simplification and clarity, we will refer to 
this ratio as our 9%/91% e-mix (for emissivity mix) in 
the rest of this article. We emphasize however that this 



is not equal to the luminosity ratio of the UV-heated to 
optically-heated component of the simulated echo spec- 
trum, which we will show in the subsequent Sections is 
of the order of ^80%/20% in the 5-38 /im range. 

We also note that this 9%/91% ratio is a direct con- 
sequence of our assumptions to have a symmetric echo 
region. Rapid size variations would lead to an alteration 
of the mixing ratio, and potentially rapid variability of 
the light echo. With only one echo spectrum, we are 
not able to address these questions, which will require 
the monitoring of one given echo region over a relatively 
short period of time. 

Ultimately, in Eq. [T^l Mdust is the only undefined pa- 
rameter, and we use it as the fitting parameter between 
the recorded spectrum and our simulated one (i.e. we use 
it to move the simulated spectrum vertically to match the 
intensity of the recorded echo spectrum) . 

6.2. Simulation results 

The resulting simulated spectrum (labelled SIMl), us- 
ing standard ISM dust (descri bed in Sec. |3]) and a 
9%/91% e-mix (described in Sec. l6.ip . is shown in Fig. [9] 
(top panel). Fitting the intensity of the simulated spec- 
trum to that of the recorded one (thick line) using 
Eq. [13] leads to an echoing dust mass of 6.6 x 10^^ g, 
or 3.3 X 10"'* M©. The dashed and dotted lines corre- 
spond to the individual contributions from UV-heated 
and optically- heated dust, respectively. In the 5-38 /xm 
range, the UV component of the simulated echo spec- 
trum accounts for 80.0% of the total luminosity of the 
simulated echo. Quite interestingly, the optical compo- 
nent influence becomes more important than that of the 
UV one beyond 50 /zm. The presence of an optical com- 
ponent especially appears essential in order to fit the 
70/im photometric measurement. It will therefore be of 
critical importance to obtain such photometric measure- 
ments beyond 50/iin in future studies of infrared echoes 
if one intends to better characterize the optically-heated 
dust component of such echoes. 

Given the total echoing dust volume found in Eq. [3 
the resulting echo region density is pccho — 640 H cm~'^ 
or 1.2 X 10~^^ g cm~^, with a column density of 2.8 x 
10^^ cm~^. Because the volume of echoing dust calcu- 



lated in Sec. 16.11 is an upper limit of the actual volume 
of echoing dust (see Sec. [3l), the value of pecho can be 
considered as a lower limit on the actual echoing dust 
density. 

The bottom panel shows the fit residual (relative to 
the recorded echo). The general trend of the recorded 
spectrum is well reproduced above ~15 ^m, as we match 
both the recorded echo spectrum as well as the 70 /xm 
photometric point. Yet, there are some discrepancies in 
the PAHs features, where our simulated spectrum is sig- 
nificantly too strong. We especially fail, at this point, 
to account for the peculiarly small 11 /im complex. The 
discrepancies around the PAHs features is reflected in 
the difference between the recorded and simulated pho- 
tometric points shown in the bottom panel of Fig. [9l 
with a ^250% overestimate of the flux in those regions. 
In the following Section, we modify our assumptions on 
the echoing dust to improve the fit below 15 /im. 

7. DISCUSSION 
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TABLE 1 

The template UV burst and optical light curve characteristics used in our simulations. 



Burst type 


Twien [K] 


Luminosity [erg s '-] 


Absorption 


Duration [days] 


Total Energy [erg] 


Optical 
UV-B 


9.4 X 10'^ 


3.8 X 10"^ 
6.4 X 10-*'' 


No 
Yes 


25 
2.5 


8.2 X 10**" 
1.4 X 10^'' 


UV 


1.5 X 10^ 


5.9 X 10*" 


Yes 


- 


- 




Fig. 9.— Simulation SIMl. Top : 9%/91% e-mix simulation 
(thin line) fitted to the recorded IR echo (thick line, in red in the 
online version) shown in Fig. [S] and associated synthetic photo- 
metric points (black rhombuses). The white squares correspond to 
MIPS 24 ^m and 70 fiva photometric measurements. The white 
rhombuses represent the Spitzer synthetic photometric measure- 
ments for our simulation. The dashed and dotted lines correspond 
to the individual emission from UV-heated and optically heated 
dust, respectively. Bottom : fit residuals relative to the recorded 
echo spectrum. 

7.1. On the influence of the dust model 

As we base our s i mulations on the du st model of 
ILi and Drailli (|2001D : IDraine and LI ()2007D . it is clear 
that resulting dust masses and densities of the echoing 
region are a direct consequence of this model. Alter- 
nate dust models, such as the ones mention in Sec. lU 
are based on different total abundances (which can vary 
along different line-of-sights), dust characteristics, mix- 
tures and dust-to-gas ratios. Because the aim of this 
paper is not to compare the goodness of different dust 
models and their ability to reproduce our recorded echo 
spectrum, but rather to show the potential of perform- 
ing such a fit and its subsequent scientific output, we 
deliberately choose t o alter the orig i nal du st model of 
ILi and Draind (|2001| ): IDraine and Lil (|2007f ) only to ex- 
plore the potential of such modifications to improve our 
fit to the recorded echo spectrum, but not in an attempt 
to solely improve their dust model. Performing a similar 
analysis to the one presented in this article, but with a 
different dust model, has a great potential but is outside 
the scope of this paper. 

7.2. The 3-20 fim features: PAHs depletion 

PAHs can be strongly affected by UV radiation, which 
can cause photoionization, dehydrogenation or photodis- 
sociation (Tielcns 2008). If those various processes are 
not directly implemented in our simulation code, we can 
nevertheless (approximatively) take them into account 
by varying the amount of various dust components in the 
code, so as to reproduce their effect on the dust composi- 
tion. Specifically, we modify the equation responsible for 



the mix in between PAHs an d quantum-heated ca rbon 
grains as follows (see Eq. 2 in ILi and Draind (J2001D ): 

CSf (a. A) = ry(a)CpAHCr,^"(«, A) + (1 -ePAH)Cf^:(a, A) 

(15) 
where C'abs'^('^' ^) ^^ ^^^ absorption cross section of the 
carbonaceous dust, C^^^{a,X) and Cf^g(a, A) the ab- 
sorption cross sections of the PAHs and quantum-heated 
graphite gr ains respective l y, ^ pah is the mixing ratio as 
defined bv ILi and Draind (|200li ). and < r/(a) < 1 is 
a function which parametrizes the PAHs destruction ra- 
tio, and that we will define below. It should be em- 
phasized that this is not equivalent to modifying the 
grain sizes distribution function J_*kait defined by 

^ uh da •^ 

iWeingartner and Draind ()2001| ) which would then im- 
pact the carbonaceous dust as a whole. PAHs are 
much more likely to be processe d by UV radiation than 
quan tum- heated graphite grains ()Luolll993 : iHunter et al.l 
|2001| ) , and we do not see any explicit reason to alter the 
latter catego ry at first. This approa ch also agrees with 
the work of iCompiegne at alJ (|2008l ) who showed that 
the dust emission in the Horsehead Nebula photodisso- 
ciation region appears to contain a smaller ratio of PAH 
vs. quantum-heated carbon grains than usual. 

In Fig. [TOl we present the result of the simulation la- 
belled SI1VI2, with a 9%/91% e-mix, and for which we 
have taken the dust mix 77(a) as follow: 



via) 






if a < 5.5A 


0.5 


ff 5.5 < a < 


1 


ff a > 8.6A 



^.6A 



(16) 



In ot her words, with a = 1.286 x N^ (|Li and Draind 
120011 ). we assume that there are no PAHs containing less 
than ~80 carbon atoms, and a decrease of 50% of the 
amount of PAHs containing from ~80 to ^^300 carbon 
atoms, with respect to the initial ISIVI dust mix. 

As expected, the lack of small PAHs has an effect 
at small wavelengths, more precisely below ~15 ^m. 
The subsequent simulated IR spectrum has a trend that 
matches better the recorded echo spectrum. The 9 /im 
gap is better reproduced, and the 6.2 /im and 7.6 /im 
features are less intense than in SLVIl. However, this 
update of the dust mix still cannot explain the small 
strength of the PAH features in the echo spectrum, 
a point which we will discuss in Sec. 17.41 The re- 
sulting total mass of echoing dust for this simulation 
is 7.0x10^^ g, equivalent to a echo region density of 
Pccho=679 H cm~'^ or 1.3 x lO"^-^ g cm~'^, with a col- 
umn density of 3.0 x 10^^ cm~^. The UV component 
contributes 83.2% to the total luminosity of the simu- 
lated echo in the 5-38 /Ltm range. 

The chosen form of 77(a) might appear somewhat ar- 
tificial. It is important to clarify here that in our code. 
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Fig. 10. — Top : Same as Fig. |9] but for simulation SIM2 : 
9%/91% e-mix with PAHs depletion as described in Eg. 1161 Bottom 
: fit residual relative to the recorded echo spectrum (full line) and 
the SIMl residual (dash-dotted line) for comparison. PAHs deple- 
tion enables to slightly better fit the background level of the emis- 
sion at small wavelengths without impacting the fit above ~15 fira. 

the UV burst spectrum and the 77(a) function are two dif- 
ferent parameters treated separately in order to fit the 
recorded IR echo spectrum. A first justification for this 
particular 77(a) function would be to assume that the 
echoing dust content locally differs from the assumed 
standard ISM dust distribution described previously, and 
that this precise echo region should not be considered to 
be pristine and unaffected by previous interactions prior 
to being h eated by the Gas A SN light. For example, 
iMicelotta e t al. (2010a) showed that PAHs with less than 
50 carbon atoms can be destroyed if encountering shocks 
with velocities greater than 100 km s^"'^. Furthermore, 
given that the standard ISM dust content is not clearly 
defined, this justification represents a potential solution. 

A second explanation for the 77(a) function is to assume 
that the dust content is directly affected by the UV burst 
via photodissociation. In this case, the form of 77(a) is not 
independent but rather a direct consequence of the UV 
burst. In SeclHl we will discuss this possibility, and show 
that indeed, it is reasonable to assume PAH destruction 
leading to the form of 77(a) showed in Eq. [16] using a 
simple theoretical model of dust destruction. 

In our simulations, wc alter the dust composition of 
both the UV-heated and the optically-heated dust. As 
we will discuss in Sec. [SI the destruction of dust by 
the UV burst is not an equilibrium process counter bal- 
anced by dust production, but rather a instantaneous 
(on an ISM dust timescale) event. It appears rather un- 
likely that any dust formation process could replenish 
the PAHs destroyed by the UV burst before the optical 
component arrival. In that sense, it appears reasonable 
to assume that the dust composition for the optical slab 
has the same composition left-over from the UV burst. 
The study of other IR echo spectra should be, in the fu- 
ture, able to shine some more light on the behaviour of 
the dust after being processed by the UV burst of the 
SN. 

7.3. The influence of PAH ionization 

As they have an ioniza tion potential of the order of 
7 eV (jLi and Drainel[200ll) , the UV burst can ionize large 
amount s of PAHs in the ec ho region. Ionized PAHs are 
known ([Bakes et al.ll200rbl ) to have weaker 3.3 /im and 
11 /xm features compared to neutral PAHs. A stronger 



ionization fraction than in the normal ISM could be 
one explanation for the weak 11 fim complex. How- 
ever, ionized PAHs also have stronger emission from 5- 
9 fiui which would increase even more the discrepan- 
cies of the fit in this zone. Since we have adopted the 
ionization function of ILi and Draind ([2001[ ) in our dust 
model ((/)ion = '/>LDoi) see their Fig. 7), we can simu- 
late the consequences of different ionization fractions. In 
SIM3, we have assumed all of the PAHs to be ionized 
(0ion = 1): and kept all the other parameters identical to 
SIM2. The resulting fit to the recorded echo spectrum 
in shown in Fig. 1111 The resulting total mass of echoing 
dust is 7.0 X lO'^^ g, equivalent to an echo region density of 
/Occho=679 H cm~'^ or 1.3 x 10~^^ g cm"'^, with a column 
density of 3.0 x 10^^ cm'^. 
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Fig. 11.— Same as Fig. [lOl but for simulation SIMS; 9%/91% 
e-mix with PAHs depletion as described in Eq. 1161 and 100% of 
PAHs ionized. 

As expected, increasing the ionization fraction in- 
creases the 6.2-f7.7 /im to 11.3/im ratio of the PAHs 
features. As a consequence, the increased PAH ioniza- 
tion fraction improves the fit in the 11.3 fim region by 
reducing the intensity of the PAH features, but also sig- 
nificantly worsen the fit for the 6-8 fim region. It can 
also be noted than even with a 100% of PAHs being ion- 
ized, the 11.3 /im feature remains strong - too strong 
with respect to the recorded echo spectrum. Hence, if 
PAHs are in principle likely to be strongly ionized by the 
passage of the UV burst, the extent to which ionization 
infiuence the shape of echo spectrum is not evident, and 
another mechanism appears to be required to explain the 
very weak PAH features. In the remaining simulations 
of this paper, we will assume the original ionization frac- 
tion 0ion = (/f'LDoi because, from our simulations, the role 
played by PAHs ionization is not clear at this stage. 

7.4. The weak 11 fim complex: a signature of 
dehydrogenation 

The lack of any reasonably strong features in the 1 1 jjlvh 
region in the recorded IR echo spectrum is quite puzzling. 
These features are usu ally associated with C-H out-of- 
plane bending modes (I Bakes et al.l I2001al : IHonv et al.l 
I2001f ). Dehydrogenation might provide an alternative 
explanation to ionization for this peculiar 11 /im com- 
plex. 

Dehydrogenation occurs when a PAH loses some H 
atoms, in our case by absorption of a UV photon. Look- 
ing closer to the recorded spectrum, one can note, along 
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with a weak 11 /im complex, very little emission around 
the 5.27 /xm and 5.7 um features ( C-H be n d/stretch com- 
bination mode, see iBakes et al.l (|2001aD : iDraine and Lil 
([2007)). Those characteristics of the recorded echo spec- 
trum point towards a possible depletion of emission orig- 
inating from C-H bonds, and suggest a strongly dehydro- 
genated state of the PAHs. If the intensity of the 11 /um 
feature suggest a rather high dehydrogenation ratio, the 
presence of the 8.6 /.tm C-H in plane bending is an indi- 
cator that not all hydrogen atoms have disappeared. 

To take dehydrogenation of PAHs into account in our 
simulations, we define the number of hydrogen atoms in 
a dehydrogenated PAH N^j : 

NJf = (5hydro xNh (17) 

with (5hydro the dehydrogenation ratio, and Nh the num- 
ber of hydrogen atoms in the un- altered PAH molecule 
containing Nc carbon atoms. 

In Fig. [121 we show the simulated spectra SIM4 where 
the amount of H atoms per PAH is 40% (top), 20% 
(center) and 0% of the original value (i.e. (Jhydro = 
0.4; 0.2; 0.0)). All the other parameters are identical 
to the simulation SIM2 shown in Fig. [TU] For the 
three cases, the resulting total mass of echoing dust 
is 7.0x10^^ g, equivalent to an echo region density of 
Pccho=679 H cm"'^ or 1.3 x 10~^^ g cm""^, with a column 
density of 3.0 x 10^^ cm~^. Note that we have assumed 
that the 5.2 /im and 5.7 /im features are functions of the 
hydrogenation ratio H/C, thus departing on this aspect 
from Draine and Lil (|2007| ). It is not clear why they did 
not have their intensities as a function of the H/C ra- 
tio, but given our recorded IR echo spectrum, we believe 
this update to be reasonable. The UV component of the 
simulated echo spectrum accounts for 83.4%, 83.3% and 
83.3% of the total luminosity of the simulated echo for 
the 40%, 20% and 0% hydrogenation simulations. The 
respective parameters of all our simulations, presented in 
Fig. m to Fig.IH are shown in Tabled 

The 11 fim region is best fitted by the 80% dehydro- 
genation curve. Our aim in SIM4 was mainly to highlight 
the effect of dehydrogenation on the simulated echo spec- 
trum. We are here, quite clearly, pushing our simulations 
to their limits. In reality, ionization of dust might play 
along with dehydrogenation and reduce this value. It will 
likely be very interesting to study how universal these 
C-H features strengths are in various light-echo spectra, 
and if they display any spatial variations, change with 
the SN type or are function of the echo region-to-SN dis- 
tance, for example. We note here that in our simulation, 
dehydrogenation seems to be more important (and plays 
a bigger role) than PAHs destruction - a rather intuitive 
fact. One would expect PAHs first to get dehydrogenated 
by loosing their surrounding H-atoms, and only then, if 
the radiation field is strong enough, to get destroyed via 
further photodissociation events inside their more robust 
carbon skeleton. 

It can be noted that compared to the SIM2 simula- 
tion, the 6-9 /im PAH features are slightly stronger in 
the SIM3 simulation. In fact, looking closely at the resid- 
ual curves, the stronger the dehydrogenation value, the 
stronger the 6-9 /xm features get. This fact most likely 
arises as we alter and artificially remove hydrogen atoms 
from PAHs, which will re-radiate more of their internal 
energy using C-C bonds. 
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Fig. 12.— Same as Fig.[lOl but for simulations SIM4: 9%/91% 
e-mix with PAHs depletion as described in Eq. 1161 and PAHs de- 
hydrogenated at 60% (top), 80% (center) and 100% (bottom). De- 
hydrogenation might be a key factor able to explain the very weak 
11 fim PAH complex. 

At this point, we are then able to reproduce closely 
the recorded echo spectrum and photometric measure- 
ments from 5 to 70 /im, with the exception of the 5-9 /im 
region, in which our simulation has too strong PAHs fea- 
tures. We do not have direct explanation for this over- 
estimate, however, it appears quite clear that PAHs are 
strongly altered in this echo region. Especially, a more 
careful treatment of PAHs dehydrogenation than the one 
adopted here would be required to better understand the 
consequences of a strong and short-lived UV burst, and 
might potentially be able to explain the peculiar 5-9 /ini 
region. 

Finally, we note that based on the classification of 
iPeeters et al.l (|2002D . our recorded echo spectrum falls 
into the so-called AB' category, based on its 7.7 /xm com- 
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TABLE 2 

Various simulations and associated parameters. 



Name Emissivity mix 



Dust mix 



''hydro 



A^dust 

10^^ fsl 



Pccho Echo region UV/opt. contribution 

[H cm~"^] column density to the total luminosity 
10i''[Hcm-2] 



SIMl 9% UV + 91% opt. fSM 0LDO1 

SIM2 9% UV + 91% opt. ISM xrj{a) "^ </.ldoi 

SIM3 9% UV + 91% opt. ISMxr?(a) "" 100 % 

SIM4 9% UV + 91% opt. ISMxr;(a) *= </.ldoi 0.4/0.2/0.0 



6.6 


640 


2.8 


80.0% / 20.0% 


7.0 


679 


3.0 


83.2% / 16.8% 


7.0 


679 


3.0 


83.8% / 16.2% 


7.0 


679 


3.0 


83.4% / 16.6% 



''As defined in Eq. \T7\ 

^ ri{a) as defined in Eq. 1161 

plex. This is a peculiar category containing only two 
objects, a non-isolated Herbig Ae Be star and a post- 
AGB star, among the 57 sources studied, for which the 
7.7 /im complex is actually composed of two features at 
7.6 /im and 7.8 /ini with roughly equal strengths. This 
association of the unusual light-echo spectrum with the 
peculiar AB' class might eventually turn out to be of in- 
terest w hen trying to u nderstand the origin of the 7.8 fim 
feature (|Peeters et al.| [2'002'). 

8. COMPARISON WITH PREVIOUS WORK 

Using a simila r app roach to the present work, 
IDwek and Arendtl (|2008D presented simulations of sev- 
eral Cas-A light echoes in the 10-40 nm range. They es- 
pecially experimented with three different types of heat- 
ing spectra : optical light curves, UV bursts, and (so- 
called) extreme UV bursts (or EUV) . Because they used 
several different bursts characteristics, spanning a range 
of luminosities and temperatures, and managed to suc- 
cessfully reproduce their echo spectra, it is interesting 
to compare our UV-B burst based on our knowledge 
of the Gas A SN type to their models. In Fig. [131 we 
plot our UV-B burst against their 'UV18' and 'EUV18' 
curves. Both are black-body spectrum with a tempera- 
ture Tuvis = 5 X Iff* K and Teuvis = 5 x 10^ K. The 
number 18 refers to the assumed SN-echo column den- 
sity, 1.5 X 10^® cm~^. Note that compared to their Fig. 6, 
both UV18 and EUV18 curves have been scaled to ac- 
count for the different distance to the light echo (160 
ly versus 199 ly), as well as by 0.11 and 0.15 respec- 
tively to match the averaged luminosity over all their 
ec ho regions for the respect ive burst model (see Table 3 
in IDwek and Arendtl (|2008l )). 

It is interesting to note that our respective burst inten- 
sities agree very well. The total luminosities of L.^^ig = 
4.2 X 10'*'^ erg s^^ and Leuvis = 5.7 x lO'''^ erg s~* are 
slightly lower but nevertheless in good agreement with 
our UV-B luminosity of L^v^b — 6.4 x lO^** erg s~^. 
The difference is mostly a consequence of the shape 
modification of our UV-B burst. The luminosity of 
our black-body UV burst prior to the modification is 
5.9 X 10^^ erg s^-^. in per fect agreement with the work of 
IDwek and Arendtl (|2CI08f ). 

Because they did n ot use data below 13 /um, 
IDwek and Arendtl (I2008D were not able to find a model 
fitting the data significantly better than the others, with 
exception to the optical ones that proved not to be able 
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Fig. 13. — Comparison of our adopted UV burst spectrum (UV- 
B, full line, in blue in the online versio n) with the UV18 (dotted 
line) and EUV18 (dashed line) spectra of lDwek and Arendtl ||2008|) 
(see their Fig. 6). The UV18 and EUV18 curves have both been 
scaled to account for the distance difference of their associated 
echo (160 ly vs 199 ly). We have taken their luminosity to be 
0.11 X lO'^^ Lp) an d 0.15 X 10^^ L© respectively, as mentioned by 
IDwek and Arendtl (12003) (see their Table 3). 

to reproduce the light echo shape. Our simulations con- 
firm these conclusions, in that the 10-40 /im range in 
a light echo spectrum is mostly influenced by the UV 
burst heating the dust. Distinguishing between various 
UV burst characteristics only becomes possible at smaller 
wavelengths, below 10 /im, where the behaviour of the 
various PAH features, as well as the global shape of the 
spectrum contains unique signatures of the intensity and 
temperature of the UV burst. In that sense, our sim- 
ulations suggest the EUV models of Dwek and Arend^ 
((2008) to be more likely, especially as energetic photons 
are required to explain the weak 11 /im PAH complex 
with dehydrogenation. As for the optical light curve, ac- 
counting for a mere 17% in the 4-38 /xm range in our 
simulations, its influen ce really becomes critical above 
50 /tm, explaining why IDwek and ArendO (|2008|) , which 
did not use 70 /im photometry data, disregarded its in- 
fluence in their flt. 

In our various simulations, we flnd an echo region den- 
sity of 640-679 H cm~' ^, relatively comparable to the 
value of -385 H cnf^ of |Dwekand Arendtj (|2008f) . Note 
that these do not, in principle, need to agree, as different 
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light echoes in different locations have been considered 
in both works. In terms of column density, we find a 
much higher value, roughly 3.0 x 10^^ H cm~^, as com- 
pared to 5 X 10^^ H cni^^. This difference arises as we 
model our echo region as containing both a UV-slab and 
an optical slab. Taking only the UV-component of our 
echo model into account to fit the recorded spectrum, we 
would find a dust mass of 4.77 x 10^^ g. Subsequently, 
only taking the UV- slab as defining the volume of the 
echo region, our simulations would lead to a column den- 
si ty of 2. Ox 10^^ H cm~^, m ore consistent with the results 
of iDwek and Areiidtl (|2008[ ). especially their Echo 6 (see 
their Table 4). 

9. DESTROYING PAHS 

In order to improve the fit of the recorded spectrum 
(especially the zone below 15 /im), we have updated the 
dust mix in our simulations using the 77(a) factor (see 
Eq. rrSI) . The form of ?y(a) as defined in Eq. [inihas been 
plugged in as an individual parameter of the simulation. 
However, if this lack of dust is linked to the UV burst, 
then 77(a) is not independent. The strength of the burst 
will determine the amount and sizes of PAHs destroyed, 
while the duration of the burst will determine the total 
column density of the echo region in which the PAHs 
are destroyed, and the relative mix of UV-heated and 
optically-heated grains. This dependency of 7;(a) can be 
used as a consistency check for our PAHs destruction 
hypothesis. 

The fate of PAHs in the ISM, and the infiu- 
ence of photodissociation, ionization and collisional 
processes has been subject to several s tudies (e.g. 
lie Page et all |200H IMicelotta et all l2010b[ ). Here, we 
co nsider dust destruction u sing the same approach 
as iSiebenmorgen and Krugell (|201d) . Let us define 
Fburst(i^), the incoming radiation fiux heating the dust ; 
Cabs('^, Nc), the PAH absorption cross section as a func- 
tion of the frequency v and the number of C-atoms Np. 
Excited PAHs will de-excite themselves by IR emission in 
a time t~l s, unless their internal energy is higher than 
a threshold value E\\^. In this case, the PAH molecule 
will be disrupted and loose one or more atoms. For a 
PAH with Nc carbon atoms, we take the threshold en- 
ergy input within 1 s for PAH destruction, i^um, to be 



^lii 



3V] 



2 



(18) 



This relation is based on the assumption that 



(19) 

with Tdis the minimum temperature for destruction, and 
Eq = 5 eV the critical (Arrhenius) energy of the atom 
being ripped off the molecule. Let us subsequently define 
X the number of dissociated atoms in a PAH of size Nc 
receiving an energy Ei^ > E\i„^ within 1 s of time, and 
ii'kin the kinetic energy of one ripped off atom. Then, 
similarly to Eq. [19] : 



^i, 



x{Eo + E^in) = 3(A^c - x)kTdi, 



(20) 



In words, as a PAH is heated by an incoming energy 
Ein higher than the threshold value £'ii,n, it will loose 
X atoms until its internal temperature decreases down 
to Tdis- Rearranging the terms and using Eq. [191 it is 



possible to express a; as a function of the incoming energy 
E'in and the PAH size Nc : 



E, 



5[eV] 



Nc 
10 



(21) 



where i?i„ can be reached by the absorption of one single 
energetic photon, or several l ess energetic photons. We 
here strictly follow the work of'Sieben morgen and Kriigell 
(pOlO) , which we refer the reader to for more details. 

At the estimated distance of the echo region, the 
amount of time required to get hit by a photon (triangle- 
symbols line in Fig. [14]) , assuming the template UV in- 
coming radiation field described in Sect. [5] is more than 
one second for PAHs with Nc < 10^. Hence, those PAH 
can only loose atoms if they are hit by one single photon 
with an energy hiy > Enm- 

PAHs with sizes above 10^ Nc will, in average, get hit 
by several photons every second, however, taking into ac- 
count the mean energy of the absorbed photons, the total 
absorbed energy will not reach the limit value E\i^ for 
the radiation field strengths considered in this work. The 
optical light curve, as considered in this study, does not 
contain photons energetic enough to be able to destroy 
PAHs. 

For PAHs with Nc < 10^, the absorbed destroying 
energy Ein,h:^>Eum P^'' PAH per second as a function of 
the PAH size is : 

/>oo 
EinM>EiiANc) ^ / Fhurst{'^)Ci,hs{l^,Nc)dl' 

(22) 
This corresponds to a number of destroying photons 
N'y.in.M>Eii,Ti absorbed per PAH per s : 



N (N \- r ., /bu.st(^)Cabs(^,iVc) 

(23) 

which finally leads to the mean time tin,/ii/>£ii,n necessary 
to absorb one destroying photon : 



tin,hiy>Eii^iNc) 



1 



Nj^inM>Eii„,{Nc) 

and the mean energy of the destroying photons 



hv; 



in,hu>Ei 



.JNc) 



Ein,hiy>En„SNc) 
Nj^in,hu>Eu^ (Nc) 



(24) 



(25) 



In Fig. [TJ] we show the time required by a PAH to 
absorb one destructive photon, tin.h,y>_Eii„, , as a function 
of PAH size Nc (filled-dots line), along with the average 
time for a PAH to get hit by a photon of any energy 
(triangle-symbols line). We use our SN burst templates 
described in Sec. [5] and assume the distance of 199 light 
years (see Sec. [5]) from Cas A to the echo region. Under 
those assumptions, we see that it takes ^^2.5 days for a 
PAH with Nc=4:50 C-atoms to be hit by one destructive 
photon with the UV burst described in Sect. [5] shining 
on the dust particles. 

We also plot in Fig. [T3] a higher-end estimate of the 
time required for a PAHs of a certain size to be de- 
stroyed via photodissociation (empty-dots line), based 
on the conservative following approach: assuming an ini- 
tial PAH size Ncfi, we find the mean time required to get 



Probing interstellar dust with infrared echoes 



13 



10°F 



10= 



10 - 



1 0^ r 



10^ r 




200 



400 
No 



600 



800 



Fig. 14. — Mean time t required for a PAH particle to get hit 
by one destructive photon (filled-dots line) and to be completely 
destroyed (empty-dots line) as a function of the number of carbon 
atoms Nc in the PAH. The triangle-symbols line shows the time 
required to get hit by a photon of any energy, including those not 
able to rip off atoms of the PAH structure. 

hit by a destructive photon (Eq. [24)) , calculate the num- 
ber X of atoms kicked off (via Eq. [2T|) given the destruc- 
tive photon mean energy (from Eq l25|) . We then start 
again with the new, smaller, PAH with Nc^i = Nq^ — x 
atoms, and iterate the process until we reach Nc,n < 20, 
a v alue below which the PAHs are phqtoliti cally unsta- 
ble (|Li and Drainell2001l:lLe Page et al.ll2003[ ). Summing 
the time of the various steps of the loop gives us the 
mean destruction time tdostruct for a PAH of size A^c,o- 
This time can be considered as an over estimate, as a 
PAH, after having lost some atoms, might not turn into 
a stable PAHs of the new size. The C-atoms will likely 
be ripped off at random position, and thus weaken the 
PAH structure. In that sense, it may be possible that 
less impacts of destroying photons are required in order 
to achieve the d estruction of the PA H - this, however, is 
not certain (see lLe Page et alll2003D . 

As seen in Fig. [T31 the destruction time increases 
rapidly with the PAH size, and it already takes ^2.5 days 
in order to destroy PAHs with an initial size of ^^300 C- 
atoms. Let us recall here that assuming all PAHs below 
80 carbon atoms and 50% of PAHs from 80 to 300 car- 
bon atoms to be destroyed lead to the SIM2 simulation. 
These very similar values suggest, under the assumptions 
of this simple PAH destruction model, that the destruc- 
tion of PAHs by the UV burst appears plausible and 
logical. It should be noted here that we have neglected 
any dust processing mechanisms other than photodissoci- 
ation, such as recombination, hydrogenatio n or the coag- 
ulatio n of small molecules on bigger ones (jHunter et al.l 
1200 It) . We do not expect these to have an impact as im- 
portant as photodissociation in shaping the echo spec- 
trum, as the reaction timescale are expected to be sig- 
nificantly larger than the duration of the UV-burst, typ- 
ically of the order of a few years for hydrogenation and 
recombi nation (in a sta ndard ISM environment, see for 
example iLe Page et al.l [2001. 2003). 

Such a dust destruction scenario also ai jpears of spe- 
cial interest when considering the fact that IKrause et al.l 
(|2008D have found strong [C IJA9850 A carbon emission 
in their optical spectrum of an echo region around Cas 
A. The [C I] lines can be excited via two main mech- 



anisms: radiative recombination and col lisional excita- 
tion by electrons. lEscalante et al.l ()199H ) suggested the 
first mechanism to explain the [C I] emission lines from 
M 42 and NGC 2024, requiring strong radiation (10^ 
to 10^ times the ISRF) heating gas with a density of 
10^ H cm~^; a radiation field intensity not so differ- 
ent from the one considered in this work. In our dust 
model, the C/H ratio in PAHs is 16.1 ppm for N^ <100, 
26.1 ppm f or Nc <200 and 36.1 ppm for Nc <500 (see 
Table 3 in iDraine and Lil (|2007D ). Dust destruction in 
the form of r]{a) described in Eq. [15] will therefore, if 
affected PAHs are entirely disrupted, increase the ISM 
gas phase abundance of C atoms by ^24 ppm. Standard 
carbon abundances in the gas phase of low to moder- 
ate density ISM have been measured to be of the order 
of 15 ppm fe.g. ICardelh et al.l[T996t ISofia et al.l [19971 . 
while iDwekl (|1997() found abundances of the order of 50- 
100 ppm in the Cold Neutral Medium. Recently, using 
a different estimation method based on the strong tran- 
sition of CII A1334 A, ISofia et al.l (|20lH ) suggested that 
the carbon gas phase abundance might be lower by ^40% 
than previous values found using methods based on a 
weak intersystem absorption transition. In any case, our 
dust destruction scenario could increase the amount of C 
atoms in the gas phase by as much as 10%-20%, which 
could in turn potentially affect the [C I] line by the same 
fac tor. If this suggest that the [C I]]A9850 A detected 
by IKrause et al.l (|2008| ) is not an unambiguous and di- 
rect confirmation of our scenario of PAHs destruction, it 
can not be ruled out that destroyed PAHs influence its 
strength. Clearly, the study of other optical light echoes 
is required to better understand the origin and the ubiq- 
uity (or not) of the [C I]A9850 A line. 

10. CONCLUSION 

We have obtained a Spitzer IRS spectrum and MIPS 
photometric measurements of an IR echo region around 
the Cas A supernova. Using simulations of an ISM dust 
mix containing both carbon and silicate quantum-heated 
and thermal grain as well as PAHs, we can, to a large 
extent, reproduce the recorded IR echo spectrum. Specif- 
ically : 

1. We find that the infrared echo spectrum can be de- 
scribed as the thermal signature from an UV and 
an optical component. Those UV and optical com- 
ponent characteristics correspond to the values of 
the type lib SN 1993J which is considered to be 
representative of the Cas A supernova based on its 
optical spectrum. Their respective contribution in 
our best fit to the echo total luminosity are ^^83% 
and ^^17% (in the 5-38 fim range). 

2. The influence of the optically- heated dust be- 
comes more important than UV-heated dust above 
~50/zm. Especially, the presence of the optical 
component in our heating spectrum is consistent 
with the 70/xm photometric measurement of the 
recorded echo. 

3. The spectrum below 15 /Ltm is better reproduced 
when removing artificially small PAHs with less 
than 80 carbon atoms, and by removing 50% of 
PAHs with 80 to 300 carbon atoms. 
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4. The especially weak 11/xm complex of PAH features 
in the recorded echo spectrum can be reproduced 
when introducing 80% of dehydrogenation for all 
PAHs in our simulation. 



Our best fit simulation implies a lower limit for the den- 
sity of this echo region of Pccho =679 H cm~'^. The 
removal of small PAHs is consistent with a simple and 
rather conservative theoretical model of PAHs destruc- 
tion via photodissociation assuming the considered UV 
burst radiation field. The strong dehydrogenation ratio 
of PAHs, as compared to the amount of PAHs destruc- 
tion, is also consistent with this picture, where PAHs 
tend to loose hydrogen atoms first. In the future, a more 
careful modelling of PAHs destruction and dehydrogena- 
tion as compared to the one adopted here could lead to 
an even better fit to the recorded IR echo spectrum. The 
relatively weak 5-9 /im PAH complex, where our simula- 
tions differ from the recorded echo spectrum by a factor 
of '^250%, remains to be explained. 

The approach we adopted in this work, in order to 
study dust processing in the ISM, has a great poten- 
tial. Our results aim at being tested towards other echo 
regions around Gas A and other SNR, and this prelimi- 
nary study shows the wealth of potential outcomes. The 
unique nature of IR echoes enables the study of other- 
wise cold and very faint ISM structures being strongly 
processed over a very short period of time. Such ISM 



clouds, potentially being initially pristine, undisturbed, 
and located away from any altering radiation sources, 
can provide unique insights in the ISM physics. Fur- 
thermore, the large amount of echoes provides us with 
an extensive number of locations and environments to 
test the simulations and conclusions presented in this ar- 
ticle. It is currently rather difficult to obtain near IR 
spectra of echoing regions, now that the Spitzer Space 
Telescope has reached the end of its duty cycle. Fortu- 
nately, there are some exciting prospects for observations 
at those wavelengt hs, such as the James Webb Space 
Telescope (JWST) (jGardner et al.l[2006( ) or the Strato- 
spheric Observatory for Infrared Astronomy (SOFIA) 
(jBeckhn et all l2007t iGehrz et all I2009D which has seen 
its first light in May 2010. These telescopes and their in- 
struments will be ideally suited to record infrared spectra 
of echoing regions, and we expect that infrared echoes 
will provide many new and unique insights into the ISM 
dust composition and chemistry, as well as into the SN 
symmetries and characterization. 

We thank the anonymous referee for his/her construc- 
tive review. Most of this work has been carried out dur- 
ing two visits of F.V. at the Max-Planck-Institut fiir As- 
tronomic in Heidelberg. F.V. thanks Alexandra Bohm 
for her help dealing with the logistic related to those vis- 
its. This research has made use of NASA's Astrophysics 
Data System. 
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